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ABSTRACT 

We calculate one-loop R-parity-violating couplings corrections to the processes 
Ti7r and bt. We find that the corrections to the — > riv de- 

cay mode are generally about 0.1%, and can be negligible. But the corrections to 
the — !■ bt decay mode can reach a few percent for the favored parameters. 
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1. Introduction 



The minimal supersymmeytic standard model(MSSM) takes the minimal Higgs 
structure of two doublets[l], which predicts the existence of three neutral and two 
charged Higgs bosons h^,H^,A^, and H^. When the Higgs boson of the Standard 
Model(SM) has a mass below 130-140 Gev and the of the MSSM are in the de- 
coupling limit (which means that H"^ is too heavy anyway to be possibly produced), 
the lightest neutral Higgs boson may be difficult to be distinguished from the neutral 
Higgs boson of the standard model(SM). But charged Higgs bosons carry a distinctive 
signature of the Higgs sector in the MSSM. Therefore, the search for charged Higgs 
bosons is very important for probing the Higgs sector of the MSSM, and will be one 
of the prime objectives of the CERN Large Hadron CoUider(LHC). 

Current bounds on charged Higgs mass can be obtained at the Tevatron, by study- 
ing the top decay t bH'^, which already eliminates some region of parameter 
space [2], whereas the combined LEP experiments gives a low bounds approximately 
mH+ > 78.6GeV at 95%CL[3]. In the MSSM, we have mH± > 120 GeV from the 
mass bounds from LEP-H for the neutral pseudoscalar A° of the MSSM (m^o > 90.5 
GeV) [4]. 

Decays of a charged Higgs boson have been studied in the hterature[5], which have 
shown that the dominate decay modes of the charged Higgs boson for large tan f3 are 
— >■ tb and rz/^, while H"^ tb, tv^ and Wh for small tan/?. For example, for 
mH+ = 250 GeV, we have Br{E^ tb) = 0.90 and Br{H+ r+u^) = 0.06 for 
tan/3 = 5, and Br{H+ tb) = 0.64 and Br{H+ r+Vr) = 0.36 for tan/3 = 30. 
Moreover, if charged Higgs boson mass mH± is very heavy, the decay of into xf Xj 
are also important [6]. 

For all these decay channels, the one-loop Electroweak, QCD and SUSY-QGD 
corrections have been studied in detail in the previous literatures, for example see 
[7]. However, those one-loop effects were studied only in the MSSM with the discrete 
multiplicative symmetry of R-parity[8], and without R-parity, the effects of one-loop 
R-parity violating couplings on the decays of charged Higgs boson have not reported in 
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the literatures so far. In this paper, we try to fill this gap and present the calculation of 
the R-parity violating effects to the process H~ bt and H~ ti7t, which arise from 
the virtual effects of R-parity Violating couplings. The most general superpotential of 
the MSSM consistent with the SU{3) x SU{2) xU{l) symmetry and supersymmetry 
contains R-violating interactions, which are given by[9] 

Here Li(Qi) and Ei{Ui, Di) are, respectively, the left-handed lepton (quark) doublet 
and right-handed lepton (quark) singlet chiral superfields, and i/1,2 are the Higgs 
chiral superfields. The indices i,j,k denote generations and a (3 and 7 are the color 
indices, and the superscript c denotes charge conjugation. The A and A' are the cou- 
pling constants of L (lepton number)-violating interactions and A" those of B(baryon 
number)-violating interactions. The non-observation (so far) of the proton decay 
imposes very strong constraints on the product of L-violating and B-violating cou- 
plings. It is thus conventionally assumed in the phenomenological studies that only 
one type of these interactions (either L- or B-violating) exists. Some constraints on 
these R-parity violating couplings have been obtained from various analysis of their 
phenomenological implications based on experiment [10]. 

2. Calculations 

The tree-level amplitudes of the two decay modes of charged Higgs boson, as 

shown in Fig. 1(a), are given by 

(0) ietan/3m^_ 

Ml = u^PlVu (2) 

y ImwSw 

for H~ — > T^v, and 

Mi°^ = Ub{m^ tan I3Pl + m*/ tan l3PR)vt (3) 

for H~ hi, where sw = sin^vK = 1 — ^wl^'z^ Pr,l = (1 ± 75)/2- 
The above amplitudes lead to the tree-level decay width of the form 



P(0) _ E|M(°)|V/^(m^_,a^,&^) 



167rm|^_ 



(4) 



where X^lM^'^)] is the squared matrix element, which has been summed the colors 
and spins of the out going particles, \{x,y,z) = {x — y — z)"^ — and s=(l,2) 
corresponds to the decays into r^v, bt, with ai = m,-, bi = 0, and a2 = rrib, 62 = nT't, 
respectively. 

Feynman diagrams contributing to the R-parity violating corrections to H~ — > 
Ti7^,ht are shown in Fig.l(b)-(c). 

We carried out the calculation in the t'Hooft- Feynman gauge and used dimen- 
sional reduction, which preserves supersymmetry, for regularization of the ultraviolet 
divergences in the virtual loop corrections using the on-mass-shell renormalization 
scheme[ll], in which the fine-structure constant and physical masses are chosen 
to be the renormalized parameters, and finite parts of the counterterms are fixed 
by the renormalization conditions. The coupling constant g is related to the input 
parameters e, mw-i and mz via g"^ = e'^/sfy and = 1 — m^/m|-. 

The relevant renormalization constants in the calculations of the processes H~ 
Ti7^, bt are defined as 

rrifo = rrif + 5mf, {f = T,t,b) (5) 

iljfO^{l + 5ZfL)hfL + {l + SZfR)HfR, {f^T,U,t,b) (6) 

tan /3o = (1 + 5^/3) tan /3. (7) 

For dZi3, we use the on-shell fixing condition[12] 

Im{nAOzo(mio)} = 0, (8) 

where n^o^o(m^o) is the renormalized self-energy for the mixing of the pseudoscalar 
Higgs boson and the boson, then we have 

SZi3 = Im{nAozo(mio)}/(mzo sin 2/3). (9) 

Apparently, there are no R-parity violating contributions to H^o^o in our case, which 
leads to SZp — 0. 

Taking into account the R-parity violating corrections, the renormalized ampli- 
tudes for H~ — > Til-r, bt can be written as 
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where 5Mj") and SM^"'^ are the vertex corrections and the counterterms, respectively. 
The calculations of the vertex corrections from Fig. 1(b)- 1(c) result in 

SM^% = (A'333)2| Y: t^.{[-7.(m?/tan/3 + m^an/3)C^ 

+ {^iml tan (3- ^2"^?/ tan (3)Co] {ml, m%- , 0, ml^,m^, ml)}PLV^ 

2 _ > 

+ Yl GrnnUr[{l,,PL)C^{ml,mjj-,0,ml,m'l^,mlJ]vA, (11) 

m,n=l 

(^332) 1 '2{Rm2fub[-lt^{^tmb cot P + mtmbt&n P)C^ 

Idv 27r^mu, L ^^=1 

( ]^]mtmb cot /? - ^2"^t "^6 tan /3) Co] {rnl , m]j- , ml , mj^ ,ml,ml) PnVt 

Ki2K-2GmnUb[l,,C^i\{ml,m\-,ml,Q,ml^,mljPRVty (12) 

i,n=l 



+ 
+ 



m,n=l 

with 



Gmn = —mb{iJ, + ^6 tan/3) tan/3i?^2-^ni ~ ^t(^t + A*tan/3)i?^2-^f 
+mtm6(l + tan2/3)i?^2^L + {[tan/3(m^ sin 2/3 - m^an/3) - m^jRi^R^ , (13) 

where Cq, are the three-point Feynman integrals [13], A^^b are soft SUSY- breaking 
parameters, /i is the higgsino mass parameter, m^^ -^^ ^ are the stop(sbottom, sstrange) 
masses, and i?*^'') are 2x2 matrix, which are defined to transform the stop(sbottom) 
current eigenstates to the mass eigenstates. 
The counterterms can be expressed as 

(c) ig tan (3mr ( bm^ 1 1 \ , . 

dM^J^ = 1= h -oZrR + -OZ^L UrPhVv, (14) 

V2m^ \mr 2 2 / 



ml 



5Ml% = 

V2', 



mbt&n/S ( + l-SZbR + ) UbPiVu 

\mb 2 2 y 

-Fmt cot /? ( + \5ZbL + ^5^ti? ) UbPRVt, 

\ mt I I I 



(15) 



Calculating the self-energy diagrams in Fig. 2, we can get the explicit expressions of 
the renormalization constants as follows: 

m^ 327r2^ 

E + B,\{ml,m\^,mX) + {R'^,f\B, + i?o](m^, m?^, m?)} , 

m=l 



(16) 



m 



^ (A'333)' X 



167r2 



(17) 



m=l 



(''^333) X 



167r2 

E{(0'[^i + ^o](o, 



m=l 



1 



+ (0'[^i + ^o](0,m?^,m2)}, 



(18) 



167r2 

5: {(i?^2)'[i?i + i?o](m^,mL,m?) +^(i?L2)'[i?i + 5o](m^, m?^, m; 



m=l 



(19) 



(A332) X 



m=l 



(20) 



87r2 



(A332)' X 



1 



m=l ^ 



I 1-^332 J X 



87r2' 



E {iRL2nB[ + B',]{mlml^,m',) +-ARLnB[ + B',]{ml ml, ml)] , 

m=l ^ ^ 

(21) 



IGtt^ 



E (A'3'32)'(i?L)'[5i + 5o](m^ mL, ml), 



m=l 



m 
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E (A332)'(i?L)'[5; + So] (ml ml,ml). 



(22) 
(23) 



E {^'LriRLnBi + Bo](ml, ml^,ml) 



m. 



2 2 



E iK32)'{RL2)'[B[ + B',] (ml ml^,ml). 



(24) 



where Bq, Bi are the two-point Feynman integrals [13], and B[{p^, ml, ml) — dBi/dp^, 
Bq{'p'^ ,m\,m\) = dBo/dp^. 

Above we have shown the expressions of the contributions from the couphngs A333 
and A332, while the ones from the couplings A332 and A331 are similar, and can 
be obtained straightforwardly by substituting the corresponding particle masses. 

Finally, the renormalized decay width is then given by 

r, = rf + (5r(") + 5v^f^ (25) 

with 

^pW ^ A^/^(m^-_^, al hi) Rel^MW-'^MW} {h = v,c). (26) 



3. Numerical results and conclusions 

We now present some numerical results for the R-parity violating effects on the 
processes H~ — >■ bt and H~ rvr- The SM input parameters in our calculations 
were taken to be aew{mz) = 1/128.8, mw = 80.419GeV and mz = 91.1882GeV[14], 
and mt = 178.0GeV[15] and m^im^) = 4.25GeV[16]. 

In our calculation, we take the running mass mb{Q) and mt{Q) evaluated by the 
NLO formula [17]: 

rribiQ) = U6{Q,mt)U5{mt,mb)mb{mb), 

mt{Q) = UG{Q,mt)mt{mt). (27) 



The evolution factor Uf is 
Uf{Q2,Qi) = 



asm: ' 

asiQi) 



^ ^ asjQi) - as(Q2) j(f) 
47r 



^ 12 ^ 8982 - 504/ + 40/^ 

33-2/' 3(33 -2/)2 ' ^ 
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In addition, in order to improve the perturbation calculations, especially for large 
tan/3, we made the following replacement in the tree-level couplings [17]: 



Amfe = ~^^~9^ tan /37(mg^ , mg^ , Mg) + ^^^i^'M tan i5I{mi^ , mj-^ , /x) 
tan/3 ^ Af2, /x) + J(4)'^K, ^2, /^) 



(30) 



with 



1 ^ 6^ ^ 

^(«' ^' ^) = _ (^2 2) ("^^^ log ^ + ^^^^ log ^ + log (31) 

where M2 is the parameter in the chargino and neutralino matrix, and in our calcu- 
lation, we always set = 200GeV. rrig is the gluino mass, which is related to M2 by 
rUg = {as{mg)/a2)M2[l8]. 

The two-loop leading-log relations[19] of the neutral Higgs boson masses and mix- 
ing angles in the MSSM were used. For mn- the tree-level formula was used. 
Other parameters are determined as follows: 
(i) For the parameters ^ ^ and A^^h in squark mass matrices 

M| = ' \ (32) 



with 



= m| m| cos 2(3{lf' - Cg sin^ 9w), 

^RR = ^}j,D + "^g + ^1 COS 213 SIV? 9w, 



Mlr = Mrl 



At- n cot 13 (q^i) ^ 



(33) 



Aj, — ij, tan P (q — b) J 

to simplify the calculation we assumed Mq — Mfj — and At — A^, and we used 
rrit^ , , At — A^ and /i as the input parameters. We also assume mj^ ^ = "^si,2 — 
2 + 500GeV, and m^^ 2 — "^ci,2 = "^ti 2 + SOOGeV. Such assuming of the relation 
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between the squark masses is done merely for simplicity, and actually, our numerical 
results are not sensitive to the squark masses of the first and second generation. 

(ii)According to the experimental upper bound on the couplings in the R-parity 
violating inter action [10], we take the relevant R-parity violating parameters as A333 = 
A332 = A331 = 0.3, A323 = — A332 = 0.9,A3i3 = — A331 = 0.9, the remainder values 
of a' and A" are set to zero. Otherwise, the numerical results may become small 
because of the cancellation among the contributions of the involved different A' and 
a" parameters. 

Fig. 3 presents the dependence of the tree level decay widths on ttih- , where we 
have included the QCD and SUSY running effects of top and bottom quark masses. 
Prom this Figure one sees that tree level decay widths get larger with the increasing 

of TTlH- ■ 

In Fig. 4 we present the relevant R-parity violating corrections to the tree-level 
decay widths as the functions of mn- ■ In general the corrections to the rP mode are 
negligible small. In fact, the maximum of the corrections to H~ — > r^v are of order 
0.1% only. For H" — > bt, when mn- > 230GeV, the corrections can be larger than 
4%. There are many dips and peaks on the curves, arising from the threshold effects 
from the vertex corrections at the threshold point ttih- = triii + ^6 • For example, 
as shown in Fig. 4 (2), at juh- = 245.9GeV, we have uih- — rn^^ + m^^ for tan P — 40, 
and the correction to H" — > bt can get its maximal value of 10%. 

Fig. 5 show the dependence of the R-parity violating corrections on mi^. In 
general, the corrections increase with the decreasing of m^^. For example, when 
= lOOGcV, the correction to H~ — >■ r^v is 0.4% for tan/3 = 4, while the one to 
H~ bt is 4% for tan (3 = 40, and when mf^ = 300GeV, above corrections are both 
about 0.1%. 

In Fig. 6 we present the R-parity violating corrections as a function of tan /3. We 
find that the corrections are relatively larger for low and high values of tan/5, re- 
spectively, while become smaller for intermediate values of tan /3, which is due to the 
fact that there are no enhanced effects from the Yukawa couplings H~bt and H~tb 
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at medium tan (3. 

In conclusion, we have calculated the R-parity violating effects on the processes 
H~ TVr and H~ hi. These corrections arise from the virtual effects of R-parity 
violating couplings. We find that the corrections to the H~ — > rz7^ decay mode are 
generally about 0.1%, and can be neghgible. But the corrections to the H~ — > ht 
decay mode can reach a few percent in our chosen parameter space. Compared to the 
SUSY-QCD or SUSY-EW corrections, the typical values of which can be over 10% [7], 
the R-parity violating effects on the process H~ — > hi are smaller, but not negligible 
in some region of the parameter space. 
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Figure 1: Feynman diagrams contributing to R-parity violating corrections to H 
Tiyrlbt): (a) tree-level diagram; (b) — (c) are one- loop vertex corrections. 
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Figure 2: Feynman diagrams contributing to renormalization constants. 
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Figure 3: Dependence of the tree level decay widths on tuh- for (1) i7 —>■ rvr, 
assuming: = -400GeV, At ^ ^ 600GeV, and m^^ = lOOGeV; (2) bt, 
assuming: = 600GeV, At — A), — SOOGeV, and m^^ = lOOGeV. mn- runs from 
121GeV to 900GeV and 188GeV to 900GeV for (1) and (2), respectively. 
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(2) 
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Figure 4: The R-parity violating corrections as functions of mn- for (1) H~ — > r^v, 
assuming: = -400GeV, At = Ab = 600GcV, and m^-^ = lOOGeV; (2) H- bt, 
assuming: fj, — 600GeV, At — Ab — SOOGeV, and m^^ = lOOGeV. itih- runs from 
121GeV to 900GeV and 188GeV to 900GeV for (1) and (2), respectively. 
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Figure 5: The R-parity violating corrections as functions of mf^ for (1) H rvr, 
assuming: ^ = -400GeV, At ^ ^ 600GeV, and m^o = 350GeV; (2) H' bt, 
assuming: pi = 600GeV, At = Ab ^ SOOGeV, and m^o = 200GeV. 
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Figure 6: The R-parity violating corrections as functions of tan /3 for (1) if tu^, 

assuming: At = A,, = 600GeV, m^o = 350GeV, and m^^ = lOOGeV; (2) H bt, 
assuming: At = Af, = SOOGeV, tuao = 600GeV, and m^-^ = lOOGeV. 
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